arXiv:astro-ph/0512314vl 12 Dec 2005 


Astronomy & Astrophysics manuscript no. ms3011 February 5, 2008 

(DOI: will be inserted by hand later) 


The effect of a strong external radiation field on protostellar 

envelopes in Orion 

J. K. J0rgensen'’^*, D. Johnstone^’^, E. F. van Dishoeck^ and S. D. Doty^ 


* Leiden Observatory, P.O. Box 9513, NL-2300 RA Leiden, The Netherlands 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street MS42, Cambridge, MA 02138, USA 
^ National Research Council Canada, Herzberg Institute of Astrophysics, 5071 West Saanich Rd, Victoria, BC, V9E 2E7, 
Canada 

'* Department of Physics & Astronomy, University of Victoria, Victoria, BC, V8P 1 Al, Canada 
^ Department of Physics and Astronomy, Denison University, Granville, OH 43023, USA 

Received 8 March 2005 / Accepted 5 December 2005 

Abstract. We discuss the effects of an enhanced interstellar radiation field (ISRF) on the observables of protostellar cores in 
the Orion cloud region. Dust radiative transfer is used to constrain the envelope physical structure by reproducing SCUBA 
850 pm emission. Previously reported *'^CO, C'^O and H 2 CO line observations are reproduced through detailed Monte Carlo 
line radiative transfer models. It is found that the '^CO line emission is marginally optically thick and sensitive to the physical 
conditions in the outer envelope. An increased temperature in this region is needed in order to reproduce the '^CO line strengths 
and it is suggested to be caused by a strong heating from the exterior, corresponding to an ISRF in Orion 10^ times stronger 
than the “standard” ISRF. The typical temperatures in the outer envelope are higher than the desorption temperature for CO. 
The C*’0 emission is less sensitive to this increased temperature but rather traces the bulk envelope material. The data are 
only fit by a model where CO is depleted, except in the inner and outermost regions where the temperature increases above 
30-40 K. The fact that the temperatures do not drop below » 25 K in any of the envelopes whereas a significant fraction of 
CO is frozen-out suggest that the interstellar radiation field has changed through the evolution of the cores. The H 2 CO lines are 
successfully reproduced in the model of an increased ISRF with constant abundances of 3-5x10“**’. 
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1. Introduction 

Giant molecular clouds are the formation sites of massive stars 
in our Galaxy with the nearby Orion molecular clouds being the 
prime candidates for detailed studies of the earliest protostellar 
stages. An interesting difference for the studies of these cores 
compared to isolated systems is the large number of OB stars in 
the immediate vicinity: the UV radiation from these early-type 
stars ionizes their surrounding material within a few parsecs, 
and also affects the thermal balance and chemistry in inter¬ 
mediate and low-mass protostellar cores distributed over much 
larger scales. As recent studies indicate that large numbers of 
solar-type stars may be formed in these regions, it is of great 
importance to address the feedback between high- and low- 
mass star formation, in particular, for comparison with low- 
mass protostars formed in relative isolation in clouds such as 
Taurus. 

In this paper we present a radiative transfer study of three 
intermediate mass protostars in Orion from the sample of 
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Ijohnstone et alJ ll2003l) . The physical properties of their en¬ 
velopes are established from ID dust radiative transfer mod¬ 
eling, and CO isotopic and H2CO line observations are ana¬ 
lyzed using Monte Carlo line radiative transfer. In particular, 
we discuss the importance of heating the protostellar envelope 
via an enhanced external radiation field and the constraints on 
thi s heating from o ptically thick CO isotopomers. 

iLi et al.l ll2Q03l) recently estimated the gas kinetic tempera¬ 
tures for a sample of pre-stellar cores in Orion using inversion 
transitions of NH3. They found that these cores had lower tem¬ 
peratures than their surroundings, which they attributed to the 
impact of the strong interstellar radiation field (ISRF) in the 
Orion region. This is reminiscent of the situation for pre-stellar 
cores in more quiescent star forming regions where tempera¬ 
ture gradients due to a standard interstellar radiation field are 
found in radiative tra nsfer models repro d ucing submil l imeter 
continuum maps (e.g.. lEvans et al l l2()0ll) . lwilson et ^(Il999l) 
suggested that the ratios of peak temperatures of optically thick 
lines with different critical densities could be used as a diagnos¬ 
tic of temperatures. Applied to cores in Orion their results also 
suggested strong temperature gradients in the region around 
Orion BN/KL with more dense gas traced by NH3 inversion 
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lines colder (~ 24 K) than less den se gas (~ 40 K) probed by 
CO lines. In contrast Iwilson et'dl found similar temperatures 
probed by the NH 3 and CO lines for cores further south from 
the Orion BN/KL region and they argued that these cores are 
dominated by internal rather than external heating. This illus¬ 
trates that an important point for the thermal balance of proto¬ 
stellar objects is whether external heating can compete with the 
central source luminosity. 

The strong interstellar radiation field may also be re¬ 
flected in the chemistry, for example by increasing the 
photodissociation and photoionization. The strength of the 
UV field can for example be probed by the emission of 
C^ and has been applied to the Orion clouds where UV 
fields enhanc ed by factors of 10 ^ - 10 ^ have been sug¬ 
gested fe.g.. iTielens & Hollenbachl Il985albt iBurton et al.[ 


1990t IStacev et all 1993^ Mookeriea et all l2003ft . The strong 
UV field will also lead to enhanced abundances of electrons 
and thus abundance decreases of molecular ions, such as 
HCO^ and N 2 H^ that otherwise work as destroyers of common 
molecules. Support for this is found within the OMCl cloud 
core l ocated immediately be hind the Orion Nebula Cluster 
where lUngerechts et al.l(ll997l) infer abundance gradients along 
the cloud, with increasing abundances of species such as HCN, 
CH 3 OH, HC 3 N, and SO toward t he Orion BN/KL region . 

The astrochemistry study by Llohnstone et al.l 1I2OO3I) con¬ 
sidered a selection of cornerstone molecules in order to quan¬ 
tify the range of conditions for which individual molecular line 
tracers provide physical or chemical information. The mor¬ 
phological study compared a variety of l ocations along the 
Integral Shaped Filament (ISF) in Orion A llBallv et allll987t 
Johnstone & Ballv . Il999l) . chosen to represent a range of phys¬ 
ical conditions including enshrouded protostars, a bright PDR 
knot, and a shock front. The main conclusion of the paper was 
that the abundances of many molecular species correlate with 
source energetics, likely a result of the importance of tempera¬ 
ture dependent desorption in maintaining gas-phase molecules. 
A significant finding of the study was the need for a warm 
outer envelope for all of the protostellar sources, as expected 
given the proximity of numerous O and B stars. The study suf¬ 
fered, however, from a simplistic treatment of the molecular 
line abundances derived from single temperature and density 
statistical equilibrium calculations and thus was unable to fully 
consider the effect of an enhanced radiation field on the outer 
regions of the embedded protostars. 

Th is paper presents a continuation of the llohnstone et al.l 
1I2OO3I) study, performing a detailed radiative transfer analysis 
of the physical and chemical properties of a subsample of the 
regions studied in that paper. We concentrate on the embedded 
intermediate mass protostellar sources, MMS 6 , MMS9, and 
FIR4 where the assumption of constant density and tempera¬ 
ture throughout the envelope is most suspect. All three sources 
have large envelope masses, Menv > IOMq. The implied dust 
temperatures, T^, for the sources are > 15 K while the '^CO 
lines have peak brightness temperatures ^ 30 K which is dif¬ 
ficult to produce unless the cores are bathed in an enhanced 
interstellar radiation field. The properties of the cores studied 
in this pa per are summarized in T ableHlwith further details pre¬ 
sented in llohnstone et al.lil2003n . Even though the focus of this 


paper is on a few selected cores in Orion, a lot of the results 
brought up in this discussion are valid beyond these specific 
sources. 

This paper is laid out as follows: sect. |2l introduces the 
general problem, including the observational indications for 
a strong interstellar radiation field, and describes the mod¬ 
eling approach, including the possible constraints from con¬ 
tinuum and line observations. Sect. |3 then discusses possi¬ 
ble refinements and implications of the models, in particular, 
the resulting CO abundance structures and the reproduction of 
H 2 CO multi-transition observations. The evolutionary implica¬ 
tions are discussed in Sect.|3] Sect.|3concludes the paper sug¬ 
gesting possible further tests and future work. 

2. Models: constraining the ISRF 

Since the protostellar cores in Orion appear warmer than ex¬ 
pected from heating by only the low luminosity central sources, 
they require the enhanced interstellar radiation field produced 
by the nearby OB stars for additional heating. In this sec¬ 
tion we present continuum and line radiative transfer models 
which describe the thermal structure of these cores and can be 
used for comparison to the observations. First, however, we de¬ 
scribe direct evidence of the contribution from external heat¬ 
ing of these cores - in particular from '^CO isotopic lines and 
dust continuum observations obtained with the James Clerk 
Maxwell Telescope in a a; 15" beam (llohnstone et all l2003t 
LlOrgensen et aklEoP^) . 

2.1. CO 3-2 and ISRF 

Fig.^compares the '^CO 3-2 spectra toward the three sources. 
Also shown is the low -mass protostar, NGC 1333-IRAS 2A 
(ilOrgensen et al.l I 2 OO 4 I 1 which has a lower luminosity (16 Lq) 
than the Orion sources but also is located at shorter distances 
(220 pc). The peak temperatures for the three Orion sources 
are significantly higher than what is observed for NGC 1333- 
IRAS 2A. Also note that the lines are relatively Gaussian with 
less significant self-absorption. Since '^CO 3-2 is optically 
thick for the column densities representative of the studied 
cores, these spectra indicate that a relatively warm foreground 
layer must be present to prevent any self-absorption dip from 
the larger scale c old envelope. Comparin g the '^CO and C*^0 
3-2 spectra from Llohnstone etTP ll2003l) provides further evi¬ 
dence: the '^CO lines are about a factor 2 broader than those of 
C'^O, suggesting that the lines are optically thick. 

If the '^CO 3-2 line is optically thick, the peak temperature 
provides a measures of the excitation temperature in the region 
of the envelope where the optical depth reaches about unity. 
For t he high densities in the se envelopes the CO gas is in LTE 
(e.g.. l.lOrgensen et alll2002|l and the excitation temperature di¬ 
rectly measures the gas kinetic temperature. Eurthermore at 
these densities the gas and dust temperatures are expected to be 
closel y coupled ('e.g.. lCeccarelli et anil996LlDotv & Neufeldl 
I 1997 I) ' and can as a first approximation be assumed identical. 
In the following section we carry out the detailed line radia¬ 
tive transfer to calculate the full non-LTE line excitation and 
explore the relevance of the dust/gas temperature decoupling. 
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Table 1. Sample of sources in Orion. For details see l.Tohnstone et al](l2003l) . 


Source 


MMS6 

MMS9 

FIR4 

^bol 

[io] 

60 

90 

400 

Projected distance to the Trapezium“ 

[pc] 

2.9 

2.3 

1.7 

S 850f,m (peak) 

[Jy beam“‘] 

7.5 

2.5 

7.5 

/(‘'’CO 3-2)*’ 

[K km s“‘] 

40.7 

59.1 

76.2 

/(C'’0 3-2)*’ 

[K km s“‘] 

4.78 

2.88 

3.97 

/(H 2 CO 3o3 - 202 )" 

[K km s“‘] 

5.33 

3.07 

9.82 

/(H 2 CO 322 - 221 )" 

[Kkm s“‘] 

0.53 

0.28 

2.96 

/(H 2 CO 5o5 - 404)" 

[Kkm s“‘] 

2.56 

0.94 

5.24 


“Assuming a distance to Orion of 450 pc. *’Line intensity, ^ - f Tmb dv. 



Fig. 1. JCMT '^CO 3-2 spectra from the three sources in Orion 
compared to NGC 1333-IRAS 2A. 


C orrecting the peak bri ghtness temperatures 7'peak=17- 
28 K djohnstone et al.L l2003l) for the Rayleigh-Jeans approxi¬ 
mation, indicates gas temperatures of rg=24-35 K, assuming 
that the emission fills the beam. Such temperatures are, how¬ 
ever, impossible to reconcile with the central star being the 
sole source of luminosity. For a spherically symmetric, opti¬ 
cally thin envelope heated by a single central source of lumi¬ 
nosity, the dust temperature distribution is approximated by: 


Tdr) ^ 60 -) '( 

\2 X 10'' cm/ \ 


Cbol 

105 


qjl 


( 1 ) 


where q is related to the power-law index p of the dust 

opacity law, k = _ /cn tv/vn)^. through q = 2/(4 -i- /?). 

dScoville & Kwarj. 1976ilpotv & Leun j.ll994 Chandler et al 


Il998t IChandler & Richeii 2000l) . Tbis would be the tempera 

ture probed by the '^CO 3-2 lines if the emission comes from 
a shell with diameter larger than the JCMT beam which is very 
optically thick in these lines. To reach a specific temperature 
Td at this radius the required luminosity is determined from 
Eq. Q: 


Lboi =6x10^ 



( 2 ) 


Therefore, to produce peak temperatures 7’peak= 17-28 K in the 
JCMT beam at 330 GHz at the distance of Orion (460 pc), cen¬ 
tral sources of L ~ 6x lO" - 4x lO'' Lq are needed, significantly 
higher than the luminosity of the observed sources. In contrast, 
the Tpeak ~ 4 K for NGC 1333-IRAS 2A requires only 15- 
20 _k^ comparable to t he bolometric luminosity of the source 
d.Ie>rgensen et al.[E002i) . It therefore appears that the only way 
to explain the '^CO 3-2 line emission for the Orion sources is 
to introduce a strong external source of heating. 


2.2. Continuum models 


the same approach as in |J^rgensen_et^ 


To establish the physical structure of each core we adopted 

d2nn2lt. using the 
]l999h to sim- 


ID radiative transfer code Dusty (llvezic et al 
ulate 850 fj.m SCUBA images and to constrain the den¬ 
sity distrib ution of each c ore by comparing to observa¬ 
tions. As in llPrgensen et al.l the opacities for coagulated dust 
grains with thin ice mantles at a density of 10® cm ^ from 
lOssenkopf & Hennind d 19941) were adopted. The density pro¬ 
file is assumed to have a single radial power-law, n - 
rioirlropP. The power-law index, p, is constrained through 
comparison between observed and modeled brightness pro¬ 
files. TheOnon_protostars_are located in the dense ridge of the 
ISF djohnstone & Ball^ Il999h complicating the extraction of 
a unique brightness profile. Greyscale images of each protostar 
are shown in Figure |2l The hashed regions in each plot were 
not included in the circular averaging procedure used to de¬ 
fine t he brightness profiles. T he flux at 850 pm for each source 
from lj^nstone etd ] d2003l) was used to constrain the overall 
dust content in the core, that is the o ptical depth at the fidu - 
cial wavel ength (see, e . g., dis cussion in llorgensen et al.Ll200: 

d200ll) we adopt the average ISRF from 
_ dl99^ with the UV field from 

im' 


Following iFvans et all 
the solar neigh borhood of lBlacl 

IPraineldl978l) added and scale this at every wavelength with a 
constant factor ,/^/isrf (i.e.,;^fisRF = 1 for the standard ISRF). 


Fig.El illustrates the differences in the brightness profiles 
produced by varying the density profile slopes and the impact 
of the external interstellar radiation field. As shown, when in¬ 
cluding the ISRF a steeper density profile is required to fit the 
data. The best fit for MMS6 is a steep r'^-envelope, whereas 
MMS9 and FIR4 have flatter ^.envelopes. 
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Fig. 2. SCUBA 850 jjm maps of the three Orion sources (from left to right; MMS6, MMS9 and FIR4). The regions which have 
been excluded when extracting the brightness profiles are hatched out. 


Table 2. Summary of the tested models. 



ISRF 

p 

rioo 

rii 

10 * 

{TY 

K 

/("'CO)" 

K km s“* 

/(C‘'0)" 

Kkm s“* 

In Figures 

MMS6 

1 

2.0 

8.0 

37 

19 

22.3 (56) 

10.2 (2.1) 

12^: dotted 


10 

2.0 

7.2 

33 

21 

27.0 (51) 

9.95(1.9) 

12^: dashed 


100 

2.0 

6.5 

30 

26 

37.4 (46) 

10.6(1.7) 

12^: dotted-dashed 


1000 

1.5 

2.4 

5.9 

31 

58.6 (37) 

17.4(1.4) 

121): dashed 


1000-D 

1.5 

2.4 

5.9 

31 

55.4 (37) 

17.2(1.4) 



1000 

1.75 

3.9 

13 

32 

56.2 (39) 

14.2(1.5) 

121): dotted 


1000 * 

2.0 

5.6 

26 

36 

53.6 (40) 

11.1 (1.5) 

|21i,b: solid 


1000-2 

2.0 

5.6 

22 

37 

50.6 (40) 

11.4(1.5) 

121): dotted-dashed 


10000 

2.0 

4.5 

21 

54 

73.4 (28) 

10 .2(1.1) 


MMS9 

1 

1.5 

1.1 

2.7 

17 

34.6 (12) 

8.30 (0.47) 



1000 * 

1.5 

0.7 

1.7 

36 

64.3 (6.7) 

7.13(0.26) 

121;: solid 


10000 

1.5 

0.45 

1.1 

57 

70.9 (3.5) 

4.92 (0.13) 


FIR4 

1 

1.5 

2.3 

5.6 

22 

50.3 (16) 

17.5 (0.61) 



1000 

1.5 

1.8 

4.4 

34 

81.3 (11) 

16.4 (0.43) 

|21l: dotted 


5000 

1.5 

1.5 

3.7 

45 

99.2 (8.5) 

14.7 (0.32) 

|21i: solid 


10000 * 

1.5 

1.3 

3.2 

51 

105.5 (6.9) 

12.9 (0.26) 



Notes: MMS6 model marked with “-D” is with gas and dust temperature decoupling taken into account. The inner and outer radii are 50 and 
15,000 AU, respectively except for the MMS6 model marked with “-2” which has an outer radius of 30,000 AU. Models marked with 
are the best fit models used in Table Models for MMS6 have a constant turbulent broadening of 0.7 km s“*, those for MMS9 and FIR4 
1.0 km s“'. “Mass-weighted temperature: {T) = f T(r)n(r)4nP' dr/ f n(r)47Tr^ dr. ''Modeled 7 = 3-2 line intensity in 15" JCMT beam with 
optical thickness toward line and source center given in parentheses. 

Fig. El compares the dust temperature profile for MMS6 
with the JCMT beam size and the results for an optically thin 
envelope. The figure shows that a single pointing toward the 
source center will barely resolve the region where the dust tem¬ 
perature starts to increase due to the impact of the ISRF. This 
radius is therefore also very well-matched to the JCMT beam, 
which will selectively pick-up the outermost region. Since the 
envelope is optically thick, the temperature will be higher in the 
innermost region than predicted by the optically thin approxi¬ 
mation, but decrease much more rapidly outside a few hundred 
AU. The luminosity derived from Eq. El above is therefore a 
lower limit to the required luminosity, further underscoring the 


need for a strong interstellar radiation field. Finally, it can be 
seen that the temperature does drop below 30 K in the inter¬ 
mediate region at radii » 1000-6000 AU, where CO depletion 
may occur. However, in the outermost region where the temper¬ 
ature increases again, CO will desorb rapidly from dust grains 
resulting in “standard” CO abundances. 

Fig. 0 compares the temperature profiles for varying 
strengths of the interstellar radiation field and internal lumi¬ 
nosities. The figure clearly illustrates the competition between 
the interstellar radiation field and the internal source luminosity 
in defining the temperature structure - from the Lint = lOOL© 
model with a standard interstellar radiation field where only the 
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Fig. 3. Fits to the images of each of the Orion cores, a) MMS6 
models with varying ISRF (tisrf of 1, 10, 100 and 1000 with 
dotted, dashed, dotted-dashed and solid lines, respectively), b) 
MMS6 models with varying density slope (p - 1.5,1.75,2.0 
with dashed, dotted and solid lines, respectively). The dashed- 
dotted line indicates a model with a density slope of 2.0 and an 
outer radius of 30,000 AU rather than 15,000 AU. c) MMS9 
model, d) FIR4 model with varying ISRF (tisrf=1000 and 
10000 with dotted and solid lines, respectively). 

Table 3. Models results for various CO abundance structures. 



ISRF 

Tit 

[K] 

/("CO)* 

Kkm 

/(C"0)* 

K km s“’ 

MMS6 

1000 


53.6 (40) 

11.1 (1.5) 


1000" 

30 K 

40.3 (31) 

4.72(1.2) 


1000 

35 K 

26.2 (27) 

2.63(1.0) 

MMS9 

1000 


64.3 (6.7) 

7.13 (0.26) 


1000" 

30 K 

54.0 (5.6) 

4.56 (0.21) 


1000 

35 K 

30.2 (4.3) 

1.88 (0.16) 

FIR4 

1000 


81.3 (11) 

16.4 (0.43) 


1000 

30 K 

71.0 (9.3) 

10.6 (0.35) 


1000 

35 K 

48.1 (7.4) 

5.68 (0.28) 


5000" 

40 K 

66.7 (6.0) 

5.83 (0.23) 


10000 

45 K 

67.0 (4.8) 

4.92 (0.18) 


For all models, a depleted abundance Xd has been used which is a 
factor 100 lower than the canonical CO abundance Xq = 10“"*. "Best 
fit model. *See footnote b. Table 2. 


outer region (r > 10,000 AU) of the envelope is affected by the 
external heating to the Lint - IT 0 model with;^fisRF = 1000 
where the interstellar radiation field dominates the temperature 
in to radii of ~ 100 AU. For such a strong ISRF even a core 
without a central source or heating (or of very low luminos¬ 
ity) remains at high temperatures > 20 K all the way through 
to its center. Currently, the fluxes of these sources shortward 
of 100 pm are not known and model fits have only been made 
to the long-wavelength tail of the SED (see Llohnstone et all 
I 2 OO 3 I and references therein for data). It is therefore not pos¬ 
sible to separate the internal luminosity of the cores from the 



Fig. 4. Temperature profiles for MMS6; models for the en¬ 
hanced (xlOOO) and standard ISRF are indicated by the solid 
and dashed lines, respectively. The dotted line indicates the 
temperature profile for an optically thin envelope from Eq. o 
corresponding Xo q - 0.4 (or /c oc v® with /? = 1). The size of 
the JCMT beam is indicated by the vertical dotted line, the CO 
evaporation temperature (see discussion in Sect. ED with the 
horizontal dotted line. 



Fig. 5. Dust temperature profiles as function of radius (lower X- 
axis) and density (upper X-axis) for a n oc ^-envelope with 
inner and outer radii of 50 and 15,000 AU, respectively and 
no=2.7xlO^ cm“^. Models are shown for luminosities of the 
inner source of 1, 10, and 100 L© (red, green and blue lines) 
and ;t'iSRF of 1, 10, 100 and 1000 (solid, dotted, dashed and 
dashed-dotted lines) as indicated in the left and right part of the 
figure, respectively. 

contribution from the external radiation field, and the tempera¬ 
ture structure at the smallest radii therefore remains poorly con¬ 
strained. Sensitive infrared observations, e.g., from the Spitzer 
Space Telescope, are needed to separate the contributions from 
the internal luminosities and the dust heated by the external ra¬ 
diation field to the total radiation emitted by the cores - such as 
previously done for giant molecular c loud cores near t he galac¬ 
tic center based on ISO-LWS data bv luis et alJfeoOll) . 
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Another problem is the determination of the outer radius of 
the cores. These cores are not formed in isolation but located 
in the Orion ridge, which is also picked up at larger scales of 
the brightness profiles as seen in Fig. 13 The ridge is found 
to have a density profile » oc on scales up to 60-300" 
djohnstone & Balm Il999l) . It would therefore seem more re¬ 
alistic to adopt an attenuated profile at the outer edge of 15,000 
AU for the ISRF su ch as also suggested b y the modeling by 
LShirlev et alJ ll2005ll and lFvans et al.l il200-5li . For a density pro¬ 
file n cc r ^ the column density from the outer edge, ry, to 
infinity is A = n(ry)ry oc r~^ ry, i.e., declining as r“'. For 
MMS6, for example, this implies a maximum extinction be¬ 
tween the outer edge at 15,000 AU and the outside radiation 
field of Ay a; 6. If the outer radius of the cores is increased the 
resulting line intensities are unchanged, however (see TableH): 
for the envelopes with density profile power-law slopes, p, of 
1.5-2.0, the main contribution to the beam averaged column 
density is located at radii corresponding to the beam size or 
smaller. The '^CO line becomes optically thick at roughly the 
same radius (and temperature) in the envelope. 

2.3. Lines 

Adopting the physical structure of each envelope from the dust 
models, the line emissi on is then modeled using the Monte 
Carlo radiative code of ISchoier et alJ 1 I 2 OO 2 I) . This code was 
be nchmarked together with o ther line radiative transfer codes 
by Ivan Zadelhoff et a ll ll2002h . The t nolecu lar data were taken 
from the database of ISchoieret alJ (l2005l) . For each model a 
turbulent broadening of 0.7-1.0 km s * was adopted, reproduc¬ 
ing the observed line widths. For the optically thick '^CO tran¬ 
sitions inclusion of a more detailed velocity field may change 
the line opacity. On the other hand, if these objects indeed are 
young and the inside-out collapse model can be used to de¬ 
scribe the velocity field or if the infall can be described by a 
power-law velocity field, the outermost regions should be close 
to static and the '^CO 3-2 lin p dominated by turbul ent broad¬ 
ening (see, e.g., discussion in LlOrgensen et al.[l2004l) . 

To first order the gas and dust temperatures have been as¬ 
sumed to be perfectly coupled. This is justified by the high 
densities - generally > 10^ cm“^ for radii < 7500 AU 
- inferred throughout most of the envelope, which leads to 
strong thermal coupling betw een the gas and dust via colli¬ 
sions ( iDotv & Neufeldi I 1997 I) . In a strong external ISRF the 
dust temperature may be less closely coupled to the gas tem¬ 
perature due to photoelectric heating and UV pumping of H 2 . 
We tested the effect o f the ISRF by calculating the gas temper¬ 
ature as described bv IPotv & Neufeldi (Il997l) . The results are 
shown in Fig. |3 As can be seen, even in these strong exter¬ 
nal radiation fields, the densities studied here lead to dominant 
collisional coupling, and thus only marginal (< a few K) dif¬ 
ferences between the dust and gas temperatures. 

For the first radiative transfer iteration, the abundances are 
assumed to be constant adopting a canonical CO abundance 
Ao = 1 X 10 and and ratios of 70 and 1950, 

respectively. The exact CO abundance will affect the deriva¬ 
tion of the strength of the ISRF somewhat as it will regulate 



Fig. 6. Comparison between dust and gas temperature for a 
standard (left) and lOOOx enhanced (right) interstellar radiation 
field. In both panels the solid line is the dust temperature. The 
gas temperature is indicated with the dashed line (assuming a 
turbulent broadening 1.6 km s“') and the dotted line (assuming 
a turbulent broadening of 0.7 km s '). 


Table 4. Best fit parameters. 


Source 

MMS6 

MMS9 

FIR4 

Tisrf 

1x10^ 

1x10^' 

1x10'* 

CO depletion, Xg/Xy, 

100 

100 

100 

Desorption temp., Tev 

30 K 

30 K 

45 K 

[P-H 2 CO] 

2.8x10^'° 

3.6x10-'“ 

5.5x10-'° 


the location in the envelope where the observed '^CO emis¬ 
sion becomes optically thick. On the other hand, as the *^CO 
3-2 emission does become optically thick the exact abundance 
structure in the interior of the envelope where depletion may 
occur does not affect the constraints on the ISRF (see also dis¬ 
cussion below). 

As shown in Fig.0 an interstellar radiation field 1000 times 
stronger than the standard field is needed to account for the 
peak temperatures of the observed *^CO spectra and fill out 
the self-absorption dip. It gives relatively Gaussian profiles and 
line widths of 2-2.5 km s * (FWHM) in agreement with obser¬ 
vations. It is worth emphasizing here the differences between 
the observations of the '^CO and C'^O lines: the C'^O lines are 
signi ficantly narrower (1.0-1 .5 km s"'; FWHM) than the '^CO 
lines ([.fohnston^ et al.Ll2003l) (seeFig.rn. This simply reflects 
the fact that the C'^O lines are optically thin and is confirmed 
by the models. Table|2]and|3summarize the results of the mod¬ 
els for each source and Table|4]lists the parameters for the best 
fit model for each source. 

3. Discussion 

In the previous section, the contribution from the external heat¬ 
ing of the cores through the interstellar radiation field has been 
determined using the established framework of dust and gas ra¬ 
diative transfer. However, some simplifying assumptions were 
made and this section discusses how these models can be re¬ 
fined by taking the importance of chemistry into account. 

3.1. CO depletion and desorption temperature 

LlOrgensen et al.lll2002l) derived the CO abundances for a sample 
of ~20 pre- and protostellar objects, among them NGC 1333- 
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Fig. 7. Modeled and observed line profiles for '^CO (left) and 
C'^O (right) 7 = 3 - 2 for MMS6 (upper), MMS9 (middle) and 
FIR4 (lower). For MMS6 '^CO models are shown for 1, 10, 
100, 10^ and 10"^ times the standard ISRF with black dashed- 
dotted, dotted, dashed, solid and dashed-3xdotted lines, re¬ 
spectively. For the remaining sources the best fit models from 
Table|3and|4]are shown. 


IRAS 2A. Clear signs of CO freeze-out were observed, with 
the highest degree of depletion found in the objects with the 
most massive envelopes. Models of the C'^O lines suggest that 
freeze-out also occurs in the intermediate mass envelopes stud¬ 
ied in this paper; adopting a constant canonical CO abundance 
of 10“^ as was done for ^^CO in the first iteration above, over¬ 
produces the observed C'^O 3-2 line intensities by about a fac¬ 
tor three. It therefore seems that, at least in the region where 
the C^^O lines are sensitive, the abundances must be lower. 
Since the C'^O 3-2 line generally is optically thin through¬ 
out the studied envelopes, the line is sensitive to the chemical 
structure in the interior of the envelope. 

The freeze-out of CO - and thus its abundance - depends 
on temperature and density in the envelope and is therefore not 
constant but rather varies with radius. A CO desorption temper¬ 
ature higher than 20 K is suggested for low-mass protostellar 
cores from Monte Carlo modeling of multi-transition data and 
direct imaging through high-resolution interferometer observa¬ 
tions. These observations show that CO evaporates off the dust 
grains in the innermost envel opes wherever the temperature be¬ 
come s higher than 35^0 K djjzirgensen et alll2OO2Hj0rgensenl 
l2004t Ijprgensen et all l2005all . Detailed ch emical models of 
the low-mass protostar IRAS 16293-2422 bv IPotv et al J il2004l) 
suggest a desorption temperature of CO of 20^0 K. 


To simulate this effect, we introduce a simple piece-wise 
constant “drop” abundance structure: we fix the abundances 
in three regions of the envelope, adopting a standard [CO/H 2 ] 
abundance of Aq = 1 x 10 in the inner and outermost region of 
the envelope where CO desorbs due to the temperature increas¬ 
ing above the CO desorption temperature T > T^y and a de¬ 
pleted abundance Xo by two orders of magnitude in the region 
where the temperature drops below Tev This model then only 
has one free parameter, Tev, which can be fit by comparison to 
the '^CO and C'^O line intensities. In fact we are not sensitive 
to the temperature, 7’ev(in), at the innermost of these radii as 
the region inside this contributes only a negligible amount of 
material to the beam averaged C'^O column density. Similar 
abundance profiles have been suggested for low-mass proto¬ 
stellar envelopes; in these envelopes the densities in the exterior 
are low so that the timescale for CO freeze-out is longer than 
the “age” of the core, which results in un-depleted abundances 
there. This is less of an issue for the Orion cores, however, since 
for these sources the densities are high enough at the edge that 
CO depletion does occur in the lifetime in the core, provided 
the temperature is low enough. 

For MMS6 and MMS9 the ^^CO and C'^O line intensi¬ 
ties are well reproduced by these structures where the abun¬ 
dance decreases when (going from the inside) the temperatures 
decreases below 30 K at «il300 AU (MMS6) and «i2300 AU 
(MMS9) and the abundance increases again when the temper¬ 
ature increases above 30 K (due to the external heating) at 
radii of «;6400 AU (MMS6) and ^5300 AU (MMS9) (Table|3. 
'^CO still becomes optically thick in the outermost envelope 
and is thus not sensitive to the depletion region, whereas the 
C^^O lines require an abundance decrease by 2 orders of mag¬ 
nitude. For both sources, models with Tev = 30 K are found 
to provide good fits to the data. A very interesting point is il¬ 
lustrated by the fits to FIR4; for this source it is not possible 
to fit the ^^CO and C'^O simultaneously with^isRp = 1000; 
for Tev = 30 K the C'^O line is overproduced by a factor 2 
by the envelope models whereas the ^^CO 3-2 line is well- 
reproduced. Increasing the CO desorption temperature to 35 K 
(i.e., increasing the size of the CO depletion zone) lowers the 
C^^O 3-2 line intensity to the observed value, but due to the 
'^CO 3-2 line being only moderately thick, the temperature 
probed by this line simultaneously drops and the line inten¬ 
sity correspondingly decreases by 35%. If ;t'isRF is increased 
by a factor 10 the envelope naturally becomes warmer, never 
dropping below 35 K. However, for such a model both ^^CO 
and C'^O line intensities are well-reproduced with a desorp¬ 
tion temperature of 45 K (corresponding to radii of 1700 AU 
and 7700 AU). This suggests a difference in the environment 
of FIR4 and the two other cores, with a stronger contribution 
of the ISRF impacting FIR4. Such a difference is in fact reason¬ 
able: FIR4 is located at about half the distance from the illumi¬ 
nating Trapezium stars compared to MMS6 (Table 1), so from 
pure geometric dilution an enhancement of the ISRF by about 
a factor of 4 is to be expected. An even larger enhancement dif¬ 
ference between FIR4 and the two other sources is possible if 
the radiation field is attenuated by any intervening cloud mate¬ 
rial. 
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Fig. 8. Modeled line intensities in “jump models” for MMS6 
where the C'^O abundance is decreased to Xd from the out¬ 
side inwards to a specific radius corresponding to an extinction 
^v.depi. and a canonical abundance Xq of 5.1 x 10“^ (correspond¬ 
ing to a CO abundance of 10 within this radius. The density 
and temperature profiles are otherwise similar to the best fit 
model for MMS6. The contours are given in steps of 1 K km s * 
ranging from 1 to 11 K km s ' as labeled. The dashed line in¬ 
dicates the observed line intensity for MMS6. The upper panel 
shows the temperature at the radius of the jump at A v.depl. ■ 


Rather than the '^CO and C'^O lines probing material with 
different levels of freeze-out, isotope selective photodissocia¬ 
tion may result in changed abundances for the two species. As 
argued bv ivan Dishoeck & Blackllll988ll . the photodissociation 
is dominated by line absorption (rather than continuum absorp¬ 
tion) so that self-shielding and mutual shielding lines of more 
abundant isotopes and H 2 are important for regulating the pho¬ 
todissociation rate of the various isotopes. This can in princi¬ 
ple affect the abundance ratios of the isotopic species, e.g., by 
decreasing the less abundant species such as C'^O relative to 
'^CO compared to the “standard isotope ratios”. Fig. |3 shows 
the predicted C'^O 3-2 line intensities for models with a low 
C'^O abundance from the outside of the envelope in to a ra¬ 
dius corresponding to a specific column density (or extinction, 
Av.depi.) in the envelope simulating photodissociation of C'^O. 
Inside this radius the C'^O abundance is increased to its canon¬ 
ical value assuming the standard isotope ratios. As shown the 
C'^O abundance has to either be lower by one to two orders of 
magnitude into an Ay.depi of about 15-20 or lower by a factor of 
about 3 into an Ay.depi. ^ 50. Photodissocation will most likely 
only affect the outermost (Ay < a few) envelope and only de¬ 
crease the abundance by a few, it is therefore not sufficient to 
invoke isotope selective photodissociation to explain the differ¬ 
ences between the material probed by '^CO and C'^O. 


3.2. H 2 CO emission 


An interesting test of the derived models is how they account 
for the emission from other molecular species. In particular, ra¬ 
tios between a few key li nes of H 2 CO have been sugg ested as a 
useful temperature probe jMangum & Woott^l 19931) and may 
therefore be particular ly sensitive to the impact of the ISRF. 
TlohnstonT assumed constant densities and temper¬ 

atures for each of the protostars in this paper and found tem¬ 
peratures of <50 K for MMS6, 30 K for MMS9 and 80-90 K 
for FIR4. For comparison the mass-weighted temperatures. 


f 


<T> = I T{r)n{r)4nr drj I n{r)47Tr dr. 




(3) 


from the derived envelope models are 36 K for MMS6 and 
MMS9 and 51 K for FIR4 (see TableH). The variations of these 
temperatures between the studied sources therefore do appear 
to reflect the differences in the external heating, with FIR4 
again being significantly warmer than the remaining sources. 

The enhanced ISRF and thus more strongly varying tem¬ 
perature profile may, however, produce a discrepancy for sen¬ 
sitive temperature probes such as the ratios of the H 2 CO lines. 
To quantify this issue, simple constant abundance models were 
constrained by comparison to the observed H 2 CO line inten¬ 
sities. Table 13 lists the best fit models for models with both 
Aisrf = 1 the best fit, enhanced ISRF for each source. It 
is clear that the models with the enhanced ISRF provide sig¬ 
nificantly better fits to the observed H 2 CO line intensities. The 
abundances are found to be remarkably similar at 3-5x10^'*'. 
For comnarison Lleirgensen et al.1 (l2005bl) found constant abun¬ 
dances of 3xl0“"^xl0^® for a sample of 18 low-mass proto¬ 
stars. 

The Orion cores are also interesting in the context of their 
H 2 CO abundance structures: enhanced H 2 CO abundances have 
been suggested for low-mass protostars in regions where ices 
evaporate from the dust grains at temperatures higher than 90- 
100 K although the eviden ce still is inconclus i ve (see discus¬ 
sions in iMaret et al.l(l2004l) and llOrgensen et al.lll2005bll ). Since 
the cores in Orion have higher luminosities, the innermost 
hot (T > 90 K) region is expected to be larger, and, despite 
their greater distances, these hot cores would be le ss diluted in 
the sin gle- dish beam than the sourc es studied by iMaret et"^ 
(l2004 and l.le>rgensen et akl ll2005bl) (diameters of 1.0 - 1.5" 
for the Orion sources compared to the < 0.5" for the low-mass 
protostars). We find, however, that the observed H 2 CO lines are 
best fit by a constant abundance or a modest abundance jump 
at best. This conclusion is only based on a few lines, however, 
some of which may become optically thick and therefore more 
sensitive to the enhanced temperatures due to the external ra¬ 
diation field. Also, chemistry in the outer envelope, especially 
freeze-out, may affec t the line intensities measured for H 2 CO 
dSchoier et allboO^l . Finally we reemphasize that it is not re¬ 
quired that these sources have a luminous internal source from 
either the dust modeling or the gas modeling, which also is 
important in this aspect. Further studies of high excitation tran¬ 
sitions at high angular resolution are needed to fully address 
these issues. 

The large spre ad in H 7 CO abundance s for the low-mass 
protostellar cores dlPrgensen et all l2005bl) may be related to 
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Table 5. Summary of best-fit constant abundance H 2 CO mod¬ 
els with differing contributions by the interstellar radiation 
field. All considered transitions are from para-H 2 CO. 


Source 

ISRF 

Abundance (para-H 2 CO) 

xL 

MMS6 

1 

2.3x10-'° 

5.2 


10^ 

2.8x10-'° 

2.9 

MMS9 

1 

2.2x10-'° 

5.1 


10^ 

3.6x10-'° 

1.8 

FIR4 

1 

2.2x10-'° 

6.7 


10“ 

5.5x10-'° 

3.0 

NGC 1333-IRAS 2A 0 

3.9x10-'° 

1.3* 


“Icr upper limit. *'For four fitted lines, including th e H 2 CO 524 - 423 
transition at 363.946 GHz Ijdrgensen et allEoOSlJt . 


CO denletion: lMaret et alJ(l2004l) and LlOrgensen et alJll200.5bll 
find that the average abundances of CO and H 2 CO are corre¬ 
lated in the outer cold envelopes of low-mass protostars and 
high resolution interferometer observations indi cate that H 2 CO 
follow s a similar abundance structure as CO jSchoier et al.[ 
I2OO4I) . H 2 CO abundances for the three Orion sources were 
therefore also calculated with drop abundance structures based 
on the CO results, i.e., varying the undepleted and depleted 
H 2 CO abundances (Xq and Xu, respectively) and adopting Te„ 
found for the fits to the CO lines. Fig. |9l shows the confi¬ 
dence plots. Interestingly a drop in abundance is not indicated 
by the observed line intensities. Constant abundance or small 
increases in abundance within the CO depletion region are al¬ 
lowed by the 2cr confidence levels for the observed sources. 


4. Evolutionary implications 

The points discussed above have interesting implications re¬ 
lated to the evolution of the cores. The first concerns their 
history: as shown in Fig. |3 even for a pre-stellar core, i.e., a 
core with no internal source of heating, the temperature never 
decreases below 25 K in such an enhanced external radiation 
field considered in this paper. Under these conditions a given 
CO molecule would desorb again almost instantaneously after 
freeze-out so that only little depletion should be expected. It is, 
however, clear that CO is frozen-out in parts of the observed en¬ 
velopes, suggesting that these regions must have been at lower 
temperatures earlier in the evolution of the cores. 

An explanation may be that the strength of the UV field 
in the Orion region has changed significantly during the evo¬ 
lution of the studied cores. Plenty of circumstantial evidence 
suggests that the Trapezium cluster, which produces the ma¬ 
jority of the local UV photons, is relative ly young llO’Delll 
l200lh . For example. IPalla & Stahleij lEooili determine the age 
for 0 ' OriB to be less than 100,000 years. The most massive 
star in the Trapezium, 0* OriC, shows consi derable variabil¬ 
ity in its X-ray emission jGagne et al.L 1 19971) a likely indica¬ 
tor of extreme youth. Most, if not all of the low-mass stars in 
the Trapezium cluster are still proplyds, despite losing mass 
from their disks due to UV induced photoevaporation at rates 



Fig. 9. jif^-confidence plots for fits to the H 2 CO lines toward 
MMS6, MMS9 and FIR4. The 68.3% (’Icr’; also grey-scale), 
90% and 95.4% (’2cr’) levels are shown. The abundance struc¬ 
tures are assumed to follow those of the best fit CO models, i.e., 
with a change in abundance from the canonical value, Xo to Xu 
in the intermediate regions where the temperatures drop below 
30 K (MMS6 and MMS9) and 45 K (FIR4). The dashed lines 
indicate constant abundance models. The minimumfor the 
models for each source is given in the lower left corner of each 
panel. Note that in all cases the best fit model has a higher Xo 
than Xo, but also that constant abundances are consistent with 
the data on the 2cr level. 


of 10 ~^ to 10~® Mgyr ' (iBallv et all Il998t Llohnstone et all 
Il998l) . Assuming that each disk has a minimum Solar nebula 
mass of ss 10“^ Mq, this implies a lifetime of < 10^ yr since the 
initiation of evaporation. Finally, the Hcr-tO [III] arcs observed 
around young low-mass stars near the high mass Trapezium 
stars trace out a wind driven bubble the s ize of which i s very 
small, again sugg esting extreme youth jO’Dell et al.[ 119971 
lBallvetal.[[T9m 

If the CO molecule is frozen-out at lower temperatures it 
may stay bound to the dust grains at the higher temperatures 
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Fig. 10. The evolution of the temperature profiles (grey lines) 
for a 100 Lg core with density profile as in Fig. |3 Each line 
corresponds to increasing ;ifisRF ranging from 100 to 10 "^ in 
logarithmic steps. The grey solid lines indicate the tempera¬ 
ture profiles for;t'isRF = 100, 1000 and 10,000. The thick black 
line in each plot indicates the radii where CO was frozen out 
forTiSRF = 100. In the lower panel the evolution of the inferred 
7’ev(out) with increasing ISRF strength is shown. 


observed at present day, such as seen in low-mass protostel¬ 
lar envelopes where CO is expected to freeze-out in the cold 
{T < 15 K) stages and evaporate as the protostar starts heating 
the envelope, increasing the temperature above » 35 -40 K. 
One explanation for this according to laborator y experiments is 
that C O migrates into a porous ice when heated jCollings et al.[ 
l2003h . It is therefore also not clear whether the determined val¬ 
ues for Tev in this paper correspond to the radii where the ice 
was formed at a lower temperature in an earlier stage or rather 
the radius where the temperature is high enough so that it evap¬ 
orates at present. In fact the difference between FIR4 on one 
hand and MMS 6 and MMS9 on the other may reflect such ef¬ 
fects: Fig. cni shows the evolution of the temperature distribu¬ 
tion in a 100 Lq envelope with a similar density profile as in 
Fig.El but for varying ISRF strengths. At the earliest stages 
CO freezes out between radii where the temperature is lower 
than 21 K, but as the ISRF strength increases the same radii 
will correspond to increasing temperatures. As the temperature 
then increases above Tev, CO desorbs and the radius for the 
outer edge of the depletion zone moves inwards. 


If this interpretation is true - the cores had CO freeze-out 
in the outer regions when formed, and then had CO evaporate 
when the OB stars “turned on” - this will also have impor¬ 
tant chemical implications in relation to the grain chemistry 
and for gas-phase species regulated by the amount of CO. The 
H 2 CO and CO correlation can easily can be understood since 
the amount of CO and related C^ an d O is the main bo ttle- 
neck for gas-phase H 2 CO production dSchoier e t alll20 04h. In 
a gas-p hase chemical model using the annroach of lOotv et alJ 
(120041) . and adopting the physical conditions (inch temperature 
and density) in MMS 6 , the H 2 CO abundances follows the CO 
abundance structure after ~ 10 ^ years. 

Although the abundance structures are similar in the Orion 
cores and in the low-mass protostellar envelopes, the grain 
chemistry might be significantly different. In both cases high 
abundances of, e.g., CO are expected in the outer envelopes, 
but whereas this is caused by low densities (and thus long 
timescales for freeze-out) in the low-mass protostellar en¬ 
velopes, it is a result of direct desorption (high temperatures) 
in the Orion protostars. This has the consequence that the 
timescale where CO is available for surface chemical reactions 
e.g., hydrogenation forming CH 3 OH and the rates of these re¬ 
actions are different in the Orion cores due to their dependence 
on temperature. 

It will naturally be very interesting to test the importance 
of any photochemistry. In particular an increased CO abun¬ 
dance provides high amount of atomic carbon, which will lead 
to higher abundances of for example CN and CS in the outer 
envelopes relative to related species such as HCN and SO. 

5. Conclusions 

We have presented an analysis of submillimeter continuum and 
line observations toward three intermediate mass protostellar 
cores in the Orion molecular cloud using detailed radiative 
transfer models. The main conclusions of this study are: 

- From simple luminosity and temperature considerations 
these sources must be subject to a strong external heating 
to reproduce observed '^CO 3-2 line temperatures and pro¬ 
files. 

- Detailed radiative transfer modeling confirms this, con¬ 
straining the ISRF enhancement to 10^ -10^ times the stan¬ 
dard ISRF from the studied sources. Different strengths of 
the ISRF between MMS 6 and MMS9 in one case and FIR4 
in the other can be explained in part due to the latter being 
closer to the Trapezium stars supplying the strong UV held. 

- The C'^O 3-2 lines are optically thin and thus not sensi¬ 
tive to the same temperature enhancement. Differences be¬ 
tween the '^CO and C'^O line widths reflect this difference 
between the optically thick and thin lines, which is well re¬ 
produced by the models. To reproduce the C'^O line inten¬ 
sities, significant CO depletion must have occurred in the 
part of the envelope where the temperature is lower than 
30-45 K. A difference in the CO desorption temperatures 
derived for the sources, however, suggest an evolutionary 
difference with the OB stars (and thus UV held) evolving 
over similar timescales as the protostellar cores themselves. 
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- Multi-transition H 2 CO observations indicate high temper¬ 
atures in the cores and their line ratios can only be repro¬ 
duced in models with a strong ISRF. Constant abundances 
of 3x10 '*’-5x10 provide a good fit to the line intensi¬ 
ties. 

This work illustrates the necessity of establishing the en¬ 
vironmental impact for the evolution of pre- and protostellar 
cores - for example in regions such as Orion where the impact 
of the radiation stars have a large impact on newly formed low- 
and intermediate mass protostars. As the above discussions also 
illustrate, detailed radiative transfer models can be used to ad¬ 
dress some of these issues, in particular if the chemistry is taken 
into account. 

Additional mid-infrared observations of the cores can fur¬ 
ther constrain the relative contributions of the internal and ex¬ 
ternal radiation field to the observed spectral energy distribu¬ 
tions a s it has previousl y been done for cores near the galactic 
center dl.is et al.[ l200llt . Further continuum observations with 
(sub)millimeter interferometers such as the SMA, CARMA 
and eventually ALMA could further confirm the temperature 
structure of the envelopes at smaller scales and in particular 
whether desorption of molecules are seen toward their centers 
due to the heating by the central protostar. Likewise in an ex¬ 
treme interstellar radiation field such as that of the Orion fur¬ 
ther studies of the chemistry in the outer envelopes and inter¬ 
vening material in the ridge are interesting to fully understand 
the impact of the ionizing sources. Understanding these pro¬ 
cesses are important as the thermal balance and pressure of the 
cores controls the conditions and outcome of the collapse and 
thus establishing the properties of the newly formed stars. 
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